Anderson localization of pairs in bichromatic optical lattices 



Gabriel Dufour and Giuliano Orso 



O 
(N 

> 

o 
O 



C/2 



I 

o 



(N 



O 



X 



Laboratoire Materiaux et Phenomenes Quantiques, Universite Paris Diderot-Paris 7 and CNRS, UMR 7162, 75205 Paris Cedex 13, France 

(Dated: March 6, 2013) 

We investigate the formation of bound states made of two interacting atoms moving in a one dimensional (ID) 
quasi-periodic optical lattice. We derive the quantum phase diagram for Anderson localization of both attrac- 
tively and repulsively bound pairs. We calculate the pair binding energy and show analytically that its behavior 
as a function of the interaction strength depends crucially on the nature -extended, multifractal, localized- of the 
single-particle atomic states. Experimental implications of our results are discussed. 



Introduction. The study of Anderson localization |[ll (AL) 
of particles in disordered or quasi-disordered potentials is a 
cornerstone in condensed matter physics. The recent observa- 
tions iH-Sl of AL in clouds of ultra-cold atoms |^ in the pres- 
ence of laser speckles or quasi-periodic optical potentials has 
opened new prospects 101 to study the rich interplay of disor- 
der and interaction effects in a highly controlled way. Prelim- 
inary experimental results [8] have been obtained for weakly 
interacting ID Bose gases subject to a quasi-periodic poten- 
tial. Theoretical studies ^] have also addressed the emer- 
gence of a Bose-glass phase in these systems. Interestingly, 
AL of light waves was recently observed ifioll in similar ID 
configurations using photonic lattices and non-linear effects 
have been investigated. 

In this Letter we address the problem of two interac ting 
particles moving in a quasi-periodic ID optical lattice 111 ll - 
1311 . We derive the phase diagram for Anderson localization of 
both attractively and repulsively bound pairs lll4ll which can 
be promptly accessed in current experiments with ultra-cold 
atoms. We calculate the pair binding energy E'b and show 
analytically that when the constituent single particle states be- 
come critical, Eb exhibits an anomalous power law exponent 
[cfr Eq fm as a function of the interaction strength. We also 
introduce a simple variational ansatz yielding very accurate 
results for the binding energy in all regimes. Since pairing 
is found to dramatically modify the localization properties of 
atoms, our approach provides also an important tool to inves- 
tigate disordered Fermi superfluids lll5l - jl7ll undergoing the so 
called BCS-Bose Einstein Condensation (BEC) crossover . 

Let us start by writing the Schrodinger equation for the sin- 
gle particle states in the absence of interaction: 

- J(j>{n + 1) - J(l){n - 1) + y(n)0(n) = £0(n), (1) 

where J is the hopping rate between nearest sites, and 
V{n) ~ A cos(27r/3n + 6) is the quasi-periodic external po- 
tential. Here A is the potential strength, (^ is a uniform (ran- 
dom) phase and /? is the ratio between the wave-vectors of 
the two laser beams (see Ref . 1 1 Si 1 191] for details). Eq.([l]i 
is generally referred as the Harper itZOll or the Aubry-Andre 
model Bill . For A = 2J the (e, /3) spectral diagram is 
the well-known Hofstadter butterfly ||22|| and for irrational 
/3 the system undergoes an AL transition. In particular, for 
A < 2 J all states are extended whereas for A > 2 J all 
states are exponentially localized with a localization length 



^ = 1/ ln(A/2 J) independent of the energy. For A = 2 J all 
eigenstates are critical and exhibit multi-fractal behavior 123 1. 
More explicitly, for large but finite system sizes L, one can 
associate to each wave-function an infinite set of fractal di- 
mensions Dq which are defined from the scaling behavior of 
l<^('^)P'' L~^iii~^\ In particular the correlation di- 
mension D2 controls the behavior of the inverse participation 
ratio near the transition point. The associated energy spec- 
trum is also characterized by a set of fractal exponents, which 
can be introduced from the dependence of the bandwidth of 
a given level on the system size, e"" — e'P oc L where 
e", eP are the eigenenergies calculated with, respectively, pe- 
riodic and antiperiodic boundary conditions. For example, for 
P = (\/5 — 1)/2 (golden ratio), the single particle ground state 
has D2 0.329 Q and 7 2.374 

In the presence of on-site (Hubbard) interactions 
U = U^^^\m,m){m,m\ between the two paiticles, the 
Schrodinger equation can be written as (iJ — i?o )!'*/') —U]^)), 
where E is the energy and _ffo is the non interacting two-body 
Hamiltonian. By applying the Green's function operator 
Ge = {E — Hf))^^ both sides of the equation, we find 
lip) = GeUI'iP)- Projecting over the state \n,7i') gives 



ilj{n,n') = U n'\GE\fn, m)tp{m, m), 



(2) 



where ip{n,n') ~ {n,n'\ip) is the amplitude of the 
two-particle wave-function. The matrix elements in 
Eq.dJI can be obtained by expressing the Green's func- 
tion operator in terms of the non interacting eigenbasis, 
Ge = Y.r,si^ ^ - es)"^|0r,(/>s)(0r>s|, whcrc (j)r are 
the eigenstates of Eq.([T]) with energy written in ascending 
order ei < 62 < ... < €l', this gives 

{n,n \GE\m,m} ^ } . (3) 

^ E - Er - Ss 

r,s 

Eq.(|2]i and ^ show that for contact interactions the two- 
body wave-function can be entirely reconstructed from the di- 
agonal terms /(m) = ■^/'("^, m). By setting n' = n in Eq.(|2|, 
we end up with the following eigenvalue problem ll25ll 



(4) 



where the kernel is KE{n,m) = {n,n\GE\'m,m) and the 
eigenvalue is A = !/[/. For values of the energy below the 



two-particle non interacting spectrum {E < 2ei), the eigen- 
values A are all negative corresponding to attractively bound 
states, since U < and the pair has a finite size. In particu- 
lar the ground state energy is obtained by varying E until the 
lowest eigenvalue A equals 1 /U. When the energy is above the 
two-particle non interacting spectrum {E > 2eL), the eigen- 
values A are all positive corresponding to repulsively bound 
states, since U > 0. In particular, a wave-function f{n) de- 
scribing an attractively bound state with energy E represents 
also a repulsively bound state with energy —E provided the 
uniform phase 6 is shifted by tt. 

For numerical stability, the irrational number /3 is usually 
expressed as the limit of a continued fraction. For definite- 
ness we take /3 = (v^ — l)/2 which can be approximated as 
/3 ~ Fj^i/Fj, where Fj are Fibonacci numbers (defined by 
Fo = 0, Fi ^ 1 and Fj = Fj_i + Fj_2 for j > 2), and j 
is sufficiently large. Finite size effects are minimized by fix- 
ing the length of the chain to L = Fj and imposing periodic 
boundary conditions. Hereafter we set J = 1 and measure all 
energies in units of the tunneling rate. 

Anderson localization of pairs. Eq.Q should be inter- 
preted as an effective single-particle Schrodinger equation for 
the center-of-mass motion of the pair In the absence of the 
quasi-periodic potential (A = 0), the solutions of Eq.® are 
delocalized Bloch states f{ri) = e*'^"/\/i, since the quasi- 
momentum k of the pair is conserved. In the opposite atomic 
limit (A — > oo), the kernel in Eq.© becomes diagonal in real 
space, KE{n,m) = 5n,m/{E — 2V{n)), implying that the 
pair is localized f{n) = Sn.na, since the two particles must be 
in the same site no to bind together 

To calculate the critical strength A = Acr where 
AL occurs, we introduce the inverse participation ratio 
Qfp = X]n=i l/l*^)!^ of the pair In the localized phase ap is 
always finite whereas in the extended phase it vanishes in the 
thermodynamic limit as ap cx L^^. For U < 0, we calculate 
ap for the ground state solution of Eq.(l4]i with fixed energy E 
and increasing values of the potential strength. We then iden- 
tify Acr as the inflection point of the obtained curve ap(A) 
1 2611 . Repeating the procedure for different E gives the quan- 
tum phase diagram shown in FigH] For vanishing interactions 
the pair breaks and Ac, — > 2 as expected for single atoms. As 
the interaction strength increases (in modulus), we see that the 
localized phase extends progressively to smaller and smaller 
values of A. This result can be better understood in the strong 
couphng regime \U\ ^ A, 1 where pairs are tightly bound 
objects with large effective mass. By expanding the kernel in 
powers of up to the third order included, we find 



Ksin, m) ~ S„ 



2V{n) 4V^(n)2+4 



£2 



E^ 



(5) 



showing that to a first approximation E ^ U [see Eq.©]. For 
tightly bound pairs we expect Ad <C 1, so we can safely ne- 
glect the term in the right hand side of Eq.®. As a conse- 
quence, the motion of the pair is governed by an equation like 




FIG. 1. Quantum phase diagram for Anderson localization of pairs 
in a ID quasi-periodic lattice. The critical potential strength is plot- 
ted as a function of the interaction (energy units J = 1). The 
dashed line corresponds to the asymptotic regime of tightly bound 
pairs where AL occurs at Acr = —2/U. The phase diagram applies 
to both attractively and repulsively bound pairs. 



([Hi, but with an effective tunneling rate Jgff = —2/E ^ 1 
and a potential strength Agff = 2A. Since AL occurs at 
Aeff = 2 Jeff, we obtain Acr = —2/E ~ —2/U, as shown 
in Fig.® with the dashed line. We also notice that the phase 
boundary between localized and delocalized molecular states 
does not depend on the choice of the angle 9. Consequently, 
the phase diagram in Fig.© applies unchanged also to repul- 
sively bound pairs for [/ > 0. 

Since in a bound state the averaged distance between 
the two atoms is finite, the density profile of atoms must 
change from extended to exponentially localized as the 
phase boundary in Fig.([T| is crossed. To see this, from 
Eqs (in and ^ we reconstruct the ground state wave- 
function, normalized to j lV'(*jj)P — 1- In Fig.(|2]i we 
plot the calculated local density rii = 2 ^„J-0(i, m)p 
(left panel) and the quasi-momentum distribution 
nu = 2^,^.„^(i,m)*V'(j,m)e*'=(^'^'' (right panel) of 
the two constituent atoms for fixed U = — 3 and for in- 
creasing values of A. We emphasize that the observed 
localization of the density profile is induced by interactions, 
as non interacting atoms would remain delocaUzed for any 
A < 2. In contrast the quasi-momentum distribution of atoms 
behaves smoothly across the phase boundary, as shown in the 
right panel of Fig.(|2]i. 

Binding energy. The pair binding energy Ei, is defined in 
the usual way from E = 2ei — Eh, where ei is the ground 
state energy for a single particle [see Eq.(|IT|)]. In Figl3](main 
panel) we show our numerical results for the binding energy 
as a function of the attractive interaction U and for increasing 
values of A (solid lines). We see that in general the quasi- 
periodic potential favors the formation of molecules, because 
interaction effects are enhanced. For |C/|A ^ 1, the molecule 
is trapped near the minimum of the external potential and we 
can treat the hopping term perturbatively. To second order 
included we find = t/ - 2A + A/ ill - A(l - cos(27r/3))). 
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FIG. 2. (Color online) Crossing the phase boundary of FigjT] Local 
density rii (left panel) and quasi-momentum distribution rik [normal- 
ized to / Ukdk / {2tt) = 2] (right panel) of atoms inside molecules 
for fixed value of U — —3 and increasing values of the potential 
strength (starting from the top). At A = Aci = 0.544 the molecular 
state becomes Anderson localized. We used 9 = n/5 and L = 233 




FIG. 3. (color online) Main panel: Binding energy of molecules 
versus interaction for increasing values of the potential strength 
A = 1 (bottom), 1.5, 2, 2.5. The asymptotic behaviors for strong 
and weak interactions, obtained from Eq.(|6j] and Eq.(|7} respectively, 
are shown with dashed lines. Diamonds correspond to the predictions 
of the variational ansatz ( I12t . Upper left panel: Aa [cfr. Eq.([8]l] and 
a\ [cfr. Eq.(|9)], as a function of the potential strength. Also shown 
is the asymptotic behavior of a\ for A ^ 1. Right panel: best vari- 
ational parameter 77 as a function of the interaction strength |?7| and 
for the same values of A as in the main panel ( A = 2 is the top curve 
near U = 0). 



from which we get 



-C/ + 2A + 2£i 



[/ - A(l - cos(27r/3)) 



(6) 



For wealc interactions \\J\ <C 1, the binding energy can be 
calculated analytically starting from the free-particle ansatz 
f{n) = 01 (ri)^ in Eqs. (O and (|4|. Taking into account that 



J2n'^r{n)cj)s (n) =^Sr.s, we find 



(7) 



where a,. = J2n are overlap functions. In par- 

ticular ai = J2n 4>i{i^)'^ is the inverse participation ratio of 
the single particle ground state. 

For A < 2 all single-particles states are delocalized and the 
overlaps vanish in the thermodynamic limit as oc L~^. For 
small Eb, the leading contribution in Eq.© comes from the 
low-lying energy states, where the overlaps become uniform, 
Uj-L c:i a, a being a constant. Going to continuous variables 
J2r ~^ L j p{'^)d£ , where p{e) is the single particle density of 
states, we can write Eq.dT]) as 1 = —Ua p{e)/{Eb + 2e), 
with e = £,. — £1. Taking into account that at low energy 
p{e) ^ j4/y^, where A is a constant that depends smoothly 
on A, we find that 



Eb^^A^a'U', 



(8) 



showing that the binding energy depends quadratically on 
the interaction strength for molecules built from delocal- 
ized states. The asymptotic behavior ^ is shown with 
dashed lines in FiglS] We see that by increasing A its 
range of validity shrinks to weaker and weaker interactions. 
At the critical point, the coefficient cr in Eq.dS]) diverges as 
(7 cx (2 — A)^(^^^2^ as shown in the upper left panel of 
FigS 

On the other hand for A > 2, the single particle spectrum is 
point-like and all the states are localized. Taking into account 
that the low-lying states have vanishing overlaps ar ~ aiSr,i, 
from Eq.(|7| we find the linear in U dependence 



EbC^lUlai, 



(9) 



typical of systems with discrete energy levels. In the atomic 
limit A ^ 1, where all states are localized within few lattice 
sites, the hopping term in Eq.([lJ can be treated perturbatively, 
yielding ai ~ 1 - 4/(A(l - cos(27r/3)))^. For weaker A, 
the inverse participation ratio decreases and vanishes at the 
critical point as ai cx (A — 2)^^ , as shown in the upper left 
panel of Figl3] 

For A = 2, where multi-fractality of the single particle 
states emerges, a fit to our numerical data reveals a power law 
behavior Eb = C\U\^, with C = 0.720 and an anomalous ex- 
ponent 5 — 1.161. In particular we find numerically that the 
overlap functions and the energy differences of the single par- 
ticle ground state scale with the system size as, respectively, 
ar cx L~^'^ and e,- — ei cx L~'', with D2 and 7 defined 
above. By substituting a,. = UrL^^^ and e,. — ei = VrL^'^ 
in Eq.(|7]i, where Ur, Vr are independent of the system size 
(provided r ^ L), we find 



1 

U 



-L-^-^^giEbL-^), 



(10) 
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where g{x) = J2'^=i^r/{x + 2Vr). Since the left hand 
side of Eq.lfTOb is L— independent, g{x) must be a power law 
g{x) cx x^^^^ from which we obtain 



7 



D2 



(11) 



which is fully consistent with our numerics. Eq.(fTTli is actu- 
ally valid for any values of A. Indeed, for A < 2 it yields 
6 = 2 because in the single particle metallic phase D2 = I 
and 7 = 2, whereas for A > 2 one finds 6 = 1 since D2 = 
in the localized phase, in agreement with Eqs.® and re- 
spectively. Eq.(|7]i nevertheless predicts a wrong numerical 
prefactor Cansatz = 0.58 at the critical point. This comes from 
the fact that for any finite interaction the molecule is localized, 
as shown in the phase diagram of Fig.([T]), whereas our naive 
free-particle ansatz remains critical. 

Variational ansatz. For finite U, we consider the following 
generalized ansatz 



— 7]d{n.nQ) 



(12) 



where r] > is a variational parameter, no corresponds to 
a site where the quasi-periodic potential takes its minimum 
value and d{n,no) represents the distance between the two 
sites. Within this class of states, the lowest eigenvalue of 
Eq.Q is given by 



1 

U 



mm- 

n 



(13) 



The obtained results are shown with diamonds in the main 
panel of Fig|3]and are in very good agreement with the full 
numerical data for all explored values of U and A. In partic- 
ular we recover the correct numerical prefactor Cvai- = 0.720 
at A = 2. In the upper right panel of Figl3]we see that 77 be- 
comes finite before the phase boundary in FiglT]is reached. To 
understand this point notice that the ansatz (fTZt applies only 
close to 77,0 so that 77 should not be confused with the inverse 
localization length of the molecule, which instead concerns 
the tails of the wave-function. 

Implications for experiments. The phase diagram ([Hi can 
be explored experimentally with sufficiently dilute bosonic 
clouds of atoms by adding a second incommensurate lattice 
to the original set-up of Ref . 11411 . Notice that the density pro- 
files in Fig |2] are calculated at equilibrium whereas in Ref.[01 
the spatial distribution of atoms is measured after letting the 
particles expand along the ID bichromatic lattice. In this case 
we expect that atoms bound into pairs will stop expanding for 
A > Acr 12711 . Finally we point out that the pair binding en- 
ergy can also be measured experimentally in optical lattices 
by rf spectroscopy 11281129 1. 
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